Impurity in a Maxwellian Unforced Granular Fluid 
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We investigate velocity statistics of an impurity immersed in a uniform granular fluid. We consider 
the cooling phase, and obtain scaling solutions of the inelastic Maxwell model analytically. First, 
we analyze identical fluid-fluid and fluid-impurity collision rates. We show that light impurities 
have similar velocity statistics as the fluid background, although their temperature is generally 
different. Asymptotically, the temperature ratio increases with the impurity mass, and it diverges 
at some critical mass. Impurities heavier than this critical mass essentially scatter of a static fluid 
background. We then analyze an improved inelastic Maxwell model with collision rates that are 
proportional to the average fluid-fluid and fluid-impurity relative velocities. Here, the temperature 
ratio remains flnite, and the system is always in the light impurity phase. Nevertheless, ratios of 
sufBciently high order moments (i'Cnpurity)/('f fluid) nia^y diverge, a consequence of the multiscaling 
asymptotic behavior. 
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I. INTRODUCTION 



Granular media are typically polydisperse. For exam- 
ple, sand and grains have a broad range of particle sizes 
and shapes. Such granular mixtures exhibit size segre- 
gation, a ubiquitous collective phenomena that underlies 
diverse processes including for example production and 
transport of powders in industry, sand dunes propaga- 
tion and volcanic flows in geophysics The "Brazil 
Nut" problem where an impurity is immersed in a uni- 
form granular media is an extreme realization of a gran- 
ular mixture as it corresponds to the vanishing volume 
fraction limit of a binary mixture. While this problem 
has been extensively studied, dynamics of such impurities 
are not fully understood |p|-[Tot. 

Understanding the velocity statistics of granular mix- 
tures is a necessary step in describing multiphase gran- 
ular flows. Experimental and theoretical studies show 
that in general, different components of a granular mix- 
ture are characterized by different typical speeds, i.e., 
granular temperatures are usually distinct [pd| . However, 
analytical treatment of this case is difhcult given the cou- 
pling between the different components of the mixture. 
An additional complication arises from the large number 
of parameters, including several restitution coefhcients, 
governing the dynamics. 

In contrast, the impurity problem is more amenable 
to analytical treatment. First, the impurity is directly 
enslaved to the fluid and second, there are fewer parame- 
ters 12 1 . In this paper, we study the impurity problem in 
the framework of the Maxwell model (l^Jl^ that assumes 
that the collision rate is uniform, i.e., independent of the 
relative velocity of the colliding particles. This simplifies 
the Boltzmann collision operator and thence this model 
is widely used in kinetic theory p^-p^. The Maxwell 
model is analytically tractable even when the collisions 
are inelastic as shown in a number of recent studies of 
uniform and polydisperse granular gases [p7|"p4[. 

We obtain analytic results for the velocity distributions 



valid for arbitrary spatial dimension and collision param- 
eters. We consider two versions of the Maxwell model. 
In the first, termed the Inelastic Maxwell Model (IMM), 
the collision rates are completely independent of the rel- 
ative velocities of the colliding particles. In the second, 
termed the Improved Inelastic Maxwell Model (IIMM), 
the collision rates arc proportional to the average relative 
velocity of the colliding particles. 

In the IMM, there are two phases separated by a crit- 
ical impurity mass m<,{rp,rq), determined by the resti- 
tution coefficients Tp and characterizing fluid-fluid 
and fluid-impurity particle collisions, respectively. When 
the impurity mass m is smaller than the critical mass, 
TO < m^, different temperatures characterize the fluid 
and the impurity. The impurity velocity distribution ex- 
hibits similar characteristics as does the fluid; in partic- 
ular, the same exponent governs the algebraic decay of 
the large velocity tail. Asymptotically, the impurity to 
fluid temperature ratio diverges as the critical mass is 
approached and is infinite when m > to* . In this regime, 
the impurity essentially scatters of an ensemble of static 
fiuid particles. The governing Lorentz-Boltzmann equa- 
tion simplifies considerably. In general, moments of the 
velocity distribution exhibit multiscaling asymptotic be- 
havior and the velocity distribution consists of replicas 
of the initial conditions. On the other hand, this phase 
transition is suppressed in the IIMM. The temperature 
ratio remains finite and the system is always in the light 
impurity phase. Nevertheless, secondary transitions af- 
fecting large velocity moments remain. 

The rest of this paper is organized as follows. In Sec. II, 
we describe the IMM and present the Boltzmann equa- 
tion for the fluid velocity distribution and the Lorentz- 
Boltzmann equation for the impurity velocity distribu- 
tion. We then determine the temperature and show that 
a phase transition occurs. The light impurity phase is 
analyzed first and then the heavy impurity phase, em- 
phasizing the behavior of the velocity distribution and its 
moments. In Sec. Ill, we discuss the IIMM and show that 



1 



the system always remains in the hght impurity phase. 
We conclude with a summary. 



II. THE INELASTIC MAXWELL MODEL 

We study dynamics of a single impurity particle in a 
uniform background of identical inelastic spheres. With- 
out loss of generality, we set the mass of the fluid parti- 
cles to unity, while the mass of the impurity is denoted 
by TO. Particles interact via binary collisions that lead 
to exchange of momentum along the impact direction. 
Collisions between two fluid particles are characterized 
by the collision parameter p, and collisions between the 
impurity and any fluid particle are characterized by the 
collision parameter q. When a particle of velocity Ui col- 
lides with a fluid particle of velocity U2 its post-collision 
velocity vi is given by 



c>P(v,t) 
di 



= J dnj duiP(ui,i) J du2P{u2,t) 



(5) 



vi = ui - (1 -_p) (g • n)n, 
vi = ui - (1 - q) (g • n)n. 



(1) 
(2) 



with g = Ui — U2 the relative velocity and n the unit vec- 
tor parallel to the impact direction. The first equation 
gives the velocity of a fluid particle and the second that 
of an impurity. The collision rules are derived by 

employing momentum conservation combined with the 
fact that in an inelastic collision, the component of the 
relative velocity parallel to the impact direction is re- 
duced by a factor equal to the fluid-fluid (impurity-fluid) 
restitution coefficient (r^). The restitution coefficients 
are related to the collision parameters via 



P = 1 - 2p, 

'q = TO — (to, -|- 1) q. 



(3) 
(4) 



Since the restitution coefhcients obey < r < 1, 
the collision parameters satisfy < p < 1/2 and 
^ 9 ^ ;;rTT' energy dissipated in each colli- 

sion is equal to A{g ■ n)^, with A — 2p{l — p) and 
A = to(1 - q)[2 ~ {m+ 1)(1 - q)] for fluid-fluid and 
impurity-fluid collisions, respectively. 

We consider a collision process where random pairs of 
particles undergo inelastic collisions with a random im- 
pact direction. This inelastic Maxwell model (IMM) is 
described by a Boltzmann equation with a uniform colli- 
sion rate. The model is a straightforward generalization 
of the classical Maxwell model. Specifically, Maxwell 
showed that for elastic "Maxwell molecules" inter- 
acting via a repulsive r~^'^'^~^^ potential in d spatial di- 
mensions, the collision rate does not depend on the mag- 
nitude of the relative velocity. 

Let P{v, t) and Q{v, t) be the normalized velocity dis- 
tributions of the background and the impurity, respec- 
tively. In this collision process, no correlations develop 
and the governing Boltzmann equation and the Lorentz- 
Boltzmann equation are 



X |(5 V - ui -I- (1 -p)(g • n)n -^(v-Ui)|, 

Jdn JduiQ{ui,t) Jdu2P{u2,t) (6) 

X |(5 V - ui (1 - q)(g • n)n - S{v - Ui)| . 

The angular integration over the impact direction is nor- 
malized to unity, J dn = 1. The collision integrals di- 
rectly reflect the collision rules. 

In writing Eqs. (^~(^, fluid- fluid and impurity- fluid 
collision rates were assumed to be the same and were 
set to unity for convenience. Hence the average num- 
ber of collisions experienced by a particle equals time. 
The results detailed in this section are exact for a mean- 
fleld collision process where, irrespective of their type, 
random pairs of particles are chosen to undergo inelastic 
collisions according to (|l|)-(^). A more realistic collision 
rates model is treated in Sec. III. 

The Boltzmann and Lorentz-Boltzmann equations 
simplify in Fourier space. Thence we use the Fourier 
transforms of the velocity distribution functions. 



F{k,t) = J dve^ "" P{v,t), 
G{k,t) = J dve''"^g(v,t). 



(7) 
(8) 



The convolution structure of Eqs. (|5|)-(^) implies that 
the collision terms factorize 

d f 

-P(k,t)-fF(k,t) = y dnP[k-p,t]F[p,t], (9) 

^G(k, t) + G(k, t) = y dn G [k - q, <] F [q, t] , (10) 

with p = (1 — p) (k • n) n and q = (1 — g) (k • n) n. These 
equations reflect the momentum transfer occurring dur- 
ing collisions. They also hint that the model is analyti- 
cally tractable. For example, expansion in powers of the 
wave number shows that moments of the velocity distri- 
butions obey closed hierarchies of evolution equations. 

A. The Temperature 

We study the freely evolving case where in the ab- 
sence of energy input the system "cools" indefinitely 
[ p5| , p6[ . The fluid temperature is deflned as the aver- 
age square of the velocity: T = ^ / rfv u^P(v, t) with 
V = |v|. From the Boltzmann equation (|5|), the fluid 
temperature evolves according to J^T — —XT with 
A = 2p{l — p) J dr\.n\. The angular integration is readily 
performed, j dxin\ = l/d, by using symmetry and the 
identity n\ + . . . + n?^ — 1. Hence the fluid temperature 
satisfies 
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d 2p{l-p) 
dt d ' 



(11) 



implying that the fluid temperature decays exponentially 
with time 



T{t)=Toe 



-2p{l-p)t/d 



(12) 



mc = 1. Thus, energy equipartition occurs on a par- 
ticular surface in the three dimensional space {m,rp,rq) 
where m(l — q)'^ = 1 — — 2p(l — p). Using relations 
(^)-(^ between the collision parameters and the resti- 
tution coefficients, energy equipartition occurs when the 
impurity mass mcq is given by 



Similarly, from the Lorentz-Boltzmann equation (|^) 
one can obtain the governing equation for the impurity 
temperature <d{t) = ^ J d\'v^Q{v,t). One finds that the 
impurity temperature is coupled to the fiuid temperature 
via a simple linear rate equation 



dt d d 



(13) 



The solution to this equation is a linear combination of 
two exponentials 

e{t) = (Go-cTo) e-(i-«')*/'' + cToe-2p(i-P)*/^ (14) 



with the constant 



(i-qr 



l-g2_2p(l-p)' 



(15) 



Therefore, there are two different regimes of behav- 
ior. When 1 — q^ > 2p{l — p), the impurity tempera- 
ture is proportional to the fiuid temperature asymptoti- 
cally, — !■ c as i ^ cxD. In the complementary region 

2p{l — p) > 1 — an extreme violation of equiparti- 
tion occurs, as the ratio of the fluid temperature to the 
impurity temperature vanishes. The impurity is very en- 
ergetic compared with the fluid and it practically sees a 
static fluid. From Eq. (p^) we find that at the transition 
point q = Y^l — 2p(l — p). Employing relations (||)-(|^) 
between the restitution coefficients and the collision pa- 
rameters we obtain the critical mass 



1 - a/(i 



-rl)/2 



(16) 



The heavy impurity phase arises when the impurity is a 
bit heavier than the fluid particles: Even when the fluid- 
fluid collisions are completely inelastic (rp = 0), the crit- 
ical mass satisfies m* > 1 -I- v^- For weakly dissipative 
fluids (r 1), the critical mass diverges, cx (1 — Tp)^^ 
(see Fig. 1). 

Note now a few features of the light impurity phase. 
First, Eq. ( |l5|) generalizes the elastic fluid (p — 0) re- 
sult c = (1 - q)/{l + q) p7|p^. That result was actually 
established for a hard sphere fluid, so at least asymptoti- 
cally both the IMM and the improved model predict the 
same impurity temperature when the fluid is elastic. Fur- 
ther, the initial impurity temperature becomes irrelevant 
and the impurity is governed by the fiuid background. 
The average energies of the impurity and fluid parti- 
cles are asymptotically equal when mO/T ^ 1, or when 



2rl 



''cq 



(17) 



As expected, this mass equals unity when Vp — Tq. Cu- 



riously, 



Icq vanishes when 



(1 + rp)/2, indicating 



that for Tq > 



y^(l + r'^)/2, energy equipartition does not 

occur, as shown in Fig. 2. However, in a generic point in 
the parameter space (m,rp,rq) equipartition does break 
down [ pTp^ . This is a signature of the dissipative and 
nonequilibriuni nature of the system. 



m* 




FIG. 1. The critical mass m* versus the restitution coeffi- 
cients Tp and rq. 



In summary, when the impurity mass is smaller than 
the critical mass m < to*, the fluid temperature governs 
the impurity temperature. Otherwise, the impurity is in- 
finitely more energetic than the fluid asymptotically. As 
i — > oo one finds 



m 



c m < TO* 

OO TO > TO* 



(18) 



At the critical impurity mass, to — to*, the solution to 
Eq. (^3|) shows that the temperature ratio diverges lin- 
early with time, — > 2p{l — p)t. Interestingly, the 
dependence on the dimension is secondary as it only sets 
the overall time scale (the transformation t ^ t/d ab- 
sorbs the dimension dependence). For example, both the 
critical mass, to*, and the temperature ratio c are inde- 
pendent of d. 
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FIG. 2. The equipartition mass rricq, given by Eq. (p^, 



versus the restitution coefficients Vp and r,. 



B. The light impurity phase 

For a light impurity, m < m, , the fluid velocity statis- 
tics govern basic characteristics of the impurity statis- 
tics. We have seen this for the temperature, as the 
initial impurity temperature is irrelevant asymptotically, 
and the ratio of the two temperatures approaches a con- 
stant 



e(t) 



:. The question arises how more detailed 
characteristics such as the impurity velocity distribution 
emerge from the fluid velocity distribution. This ques- 
tion can be answered in detail in one dimension where 
the explicit scaling solution for the fluid velocity distri- 
bution is known PQ |. Below we primarily focus on the 
one-dimensional case and derive an explicit exact solu- 
tion. We then briefly comment on the higher-dimensional 
case, where the solution is very cumbersome. 

The velocity distribution of the fluid approaches a scal- 
ing form asymptotically, P{v,t) T'^/'^V {vT'^/'^). 
Following the earlier treatment of the fluid velocity distri- 
bution [|l^,|^-|22j , we employ the Fourier transform tech- 
nique. The fluid Fourier transform is thus characterized 
by the scalin g fo rm F{k,t) = f (|fc|T^/^). The govern- 
ing equation flT]] for the corresponding scaling function 
is obtained from ^ to give 

- p(l - p)zf'{z) + f{z) = f{pz)f{z - pz), (19) 



and the solution is /(z) = (1 -I- z) e^^ |20|. 

The behavior of the fluid suggests that the impu- 
rity velocity distribution also approaches a scaling form 
Qiv,t) r-i/2Q Tj^gjj ^i^g corresponding 

Fourier transform reads G{k,t) — g (|fc| T^/^). From the 
Lorentz-Boltzmann equation (p^, this scaling function 
satisfies the linear equation 

~ p(l - p)zg'{z) + g{z) = (1 + az) e""^ g{qz) (20) 

with a = \ — q. The fluid scaling function is a combina- 



tion of ; 



with n — Q and n = 1. We therefore seek a 



similar solution of Eq. (C 



E 

n=0 



AnZ^ 



(21) 



Inserting (21) into (|2C|), the exponential factors cancel. 
Equating terms of the order z", yields the following re- 
cursion relation for the coefflcients: 



'{l-q)-p{l-p) 



\ — — np{l — p) 



(22) 



From the normalization g(0) = 1 we get Aq = I. The 
next two coefficients are Ai — 1, A2 — (1 — c)/2. Using 
these values, one computes the small z expansion of the 
solution (|l]): g{z) = 1 - cz^/2 + 0{z^), and the first 
non-trivial coefficient is indeed consistent the definition 
of the Fourier transform. 

The Fourier transform can be inverted to obtain the 
impurity velocity distribution function explicitly. The 

■J-; the inverse 



inverse Fourier transform of e '^^ is — 



7r K-^+w-^ 



transforms of z"e~^ can be obtained using successive dif- 
ferentiation with respect to k. In general, the velocity 
scaling function is 



n=0 



(23) 



In other words, one can express the solution as a combi- 
nation of powers of Lorentzians 



9 . . 

Q{w) = - V B„ 



1 



1 



(24) 



The coefficients Bn are linear combinations of the coeffi- 
cients ^fc's with k < n + 1, e.g., B2 = 1 — 3^2 + 3^3 
and B-s = 4:A2 - 24 A3 -I- 6OA4. The first squared 
Lorentzian term dominates the tail of the velocity distri- 
bution Q{w) ^ V{w) ~ w~^, as w ~^ 00. This algebraic 
behavior prevails for u; 3> 1. Therefore, the impurity has 
the same extremal velocity statistics as the fluid. 

An interesting aspect of this solution is that for partic- 
ular values of the collision parameter g, the infinite sum 
terminates at a finite order. Of course, when qi = p then 



Q(w) = V{w) 



1 



7r(l-hu;2) 



2\2 ■ 



(25) 



This is reflected by the vanishing coefficient A2 ~ 0. The 
multiplicative recursion relation (|2^) then shows that 
An — for all n > 2. Furthermore, for q^l — 92) = 
p(l — p), one has A3 — 0, and the velocity distribution 
(p4|) contains only two terms 



Qiw) 



I-3A2 



4^2 



(26) 



with A2 = q/{l + 2q). Similarly, the solution can be a 
flnite sum with k terms when p and are related via 



4 



— Qk) = p{^—p) thereby imposing A„ = for n > k. 
For given fluid collision parameter p, the above relation 
has the solution qk only for sufficiently small k. Thus, 
there are a few special values of the impurity collision 
parameter qk for which the scaled impurity velocity dis- 
tribution is a linear combination of simple rational func- 
tions (l -I- ui^) " with n = 2, . . . ,k + 1. 

While the scaling functions underlying the impurity 
and the fluid are similar, more subtle features may dif- 
fer. In particular, the full time dependent behavior, as 
characterized by the moments of the impurity distribu- 
tion, exhibits rich behavior. Let i„(t) — J dvv'^ P{v^t) 
be the moments of the fluid velocity distribution. Mul- 
tiplying the equations (||)-(||) by and integrating, the 
moments obey the recursive equations 

-Ln + anLn = ^ ( " - pY'-'L^Ln-j, (27) 



J=2 
ri-2 



dt 
with 



-M„ + 6„A/„ = ^ ( . J gJ"(l - qr-=M,Ln-j, (28) 
j=o ^-'^ 



a„(p) = l-p"-(l-p)", bn{q)^l 



(29) 



Asymptotically, the fluid moments decay exponentially 
according to [pTf 



L„(i)cxe-'^"(P)*. 



(30) 



Using this asymptotics we analyze the behavior of the im- 
purity moments. The second moment, i.e. the impurity 
temperature, was already shown to behave similar to the 
fluid temperature when 02 (p) < b2{q). The fourth mo- 
ment behaves similarly to the fourth moment of the fluid 
when 04 (p) < b^^q), and generally the first n moments 
are proportional to each other, M2 oc L2, . . . , M2„ oc L2„, 
when a2k{p) < ^2/0 (9) for fc = 1, . . . , n. In particular, all 
respective impurity and fluid moments are proportional 
to each other when a2k{p) < b2k{q) is valid for every k. It 
is easy to see that above inequalities hold when p + q < 1, 
which is always obeyed when m < 1. 

However, in the complementary parameter range, 
o-nip) > bn{q) for sufficiently large n. Such moments 
are no longer governed by the fluid and 



Mn{t) OC e 



(31) 



The corresponding moment ratio diverges asymptoti- 
cally: M„/L„ — > 00 as i ^ 00. Interestingly, the same 
behavior (|3^) is found in the heavy impurity phase, as 
will be shown below. Therefore, the two phases are not 
entirely distinct. A series of transitions affecting mo- 
ments of decreasing order occurring at increasing masses. 



TOl > 7712 > 



(32) 



signals the transition to the heavy impurity phase. When 
m > rrim the ratio M2kli^ik diverges asymptotically for 



all k > n. This generalizes the second moment transition 
occurring at 7711 = m^. The transition masses 



,)2" + (l + rp)2 



l-i[(l-rp)2" + (l + rp)2"] 



(33) 



are found from q^" = p^" + (1 — p)^" and Eq. (jij). In par- 
ticular, rrioo =lim„^oo"in = (1 + rp + 2rq)/(l - rp), and 
sufficiently light impurities (m < 777oo) mimic the fluid 
completely. Additionally, since 7n.„ > 77X00 > 1, the im- 
purity must be heavier than the fluid for any transition 
to occur. 

The above analysis of the light impurity phase suggests 
that in higher dimensions, the impurity velocity distribu- 
tion might be similar to that of the fluid. In higher di- 
mensions we again assume that the Lorentz-Boltzmann 
equation admits a scaling solution. The corresponding 
equation in Fourier space ( p^ ) then considerably simpli- 
fles. Following the earlier treatments of the fluid case 
pl| , p2| we extract the high-energy tail from the small-fc 
behavior of the Fourier transform of the scaled velocity 
distribution. The outcome is that both velocity distribu- 
tions have the same algebraic high-energy tail 



Q{w) ^ V{w) ~ w 



(34) 



with the exponent a calculated in |2l|j2^. Such analy- 
sis also yields the ratio of the prefactors governing this 
algebraic decay. 



C. The heavy impurity phase 

When the mass of the impurity is equal to or larger 
than the critical mass, 777 > 777*, the velocities of the fluid 
particles are asymptotically negligible compared with the 
velocity of the impurity. Hence, in the i — > 00 limit, fluid 
particles become stationary as viewed by the impurity. 
Therefore, one can set U2 = in the collision rule (0): 



V = u — (1 — g) (u • n) n. 



(35) 



This process is somewhat analogous to a Lorentz gas . 
However, in the granular impurity system, a heavy par- 
ticle scatters of a static background of fighter particles, 
while in the Lorentz gas the scatterers are infinitely mas- 
sive. Despite this difference, the mathematical descrip- 
tions of the two problems are similar. Specifically, the 
collision rule for the inelastic Lorentz gas is ob- 

tained from ( p5| ) by a mere replacement of the factor 
{1-q) with (l + r-g). 

Let us first consider the one-dimensional case where 
an explicit solution of the velocity distribution is possi- 
ble. Setting 772 = in the delta function in the Lorentz- 
Boltzmann equation (||) , integration over the fluid veloc- 
ity U2 is trivial, J du2 P{u2,t) = 1, and integration over 
the impurity velocity 77 1 gives 
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^^Q{v,t) + Q{v,t) ^ ig t 



(36) 



This equation can be solved directly by considering the 
stochastic process the impurity particle experiences. In a 
sequence of collisions, the impurity velocity changes ac- 
cording to vq qvo ^ q^vo ^ ■ • ■ with vq the initial ve- 
locity. After n collisions the impurity velocity decreases 
exponentially, Vn = q"'vo- The collision rate is unity, 
and hence, the average number of collisions experienced 
till time t equals t. Furthermore, the collision process is 
random, and therefore, the probability that the impurity 
undergoes exactly n collisions up to time t is Poissonian 
t"e~*/"-'- Thus, the velocity distribution function reads 



Q{v,t) 



(37) 



where Qo{v) is the initial velocity distribution of the im- 
purity. Indeed, one can check that this properly normal- 
ized solution satisfies Eq. (p6|). 

Interestingly, the impurity velocity distribution func- 
tion is a time-dependent combination of "replicas" of the 
initial velocity distribution. Since the corresponding ar- 
gument is stretched, compact velocity distributions dis- 
play an infinite set of singularities, a generic feature of 
the Maxwell model 

The impurity velocity distribution exhibits interesting 
asymptotic behaviors. Consider for simplicity the uni- 
form initial velocity distribution: Qoiv) — 1 for \v\ < 1/2 
and Qo{v) — otherwise. The solution (|3^) reduces to 
a finite sum, with n < N = ^^^j^^- In the physically 
interesting limits t — s- oo and w — > 0, the sum on the 
right-hand side of Eq. ( |37| ) simplifies to respectively e*/"? 
or (t/q)^ /N\ when the number of terms is above or below 
the threshold value N — t/q. The magnitude of Q{v,t) 
above the threshold greatly exceeds the magnitude be- 
low the threshold, so Q{v,t) appears to approach a step 
function. A refined analysis shows that the width of the 
front widens diffusively so the front remains smooth al- 
though its relative width vanishes. Specifically, one finds 
the following traveling-wave like scaling solution 



Q{v,t) 



with the following wave form and coordinate 



\/7r 



dx ( 



(38) 



(39) 



Note, however, that the large velocity tail {ij oo), ig- 
nored in Eq. ( |39|) , provides actually the dominant contri- 
bution to the moments. This is an unusual traveling wave 
form in the sense that the argument is the logarithm of 
the velocity rather then the velocity itself, a refiection of 
the exponentially decaying velocity. 

In contrast to the velocity distribution, the moments 
Mn{t) ~ J dvv"'Q{v, t) exhibit a much simpler behavior. 



Indeed, from Eq. ( |36| ) one finds that every moment is 
coupled only to itself, -^Afn _=_-(l - g")M„. Solving 



this equation we recover Eq. (31); in the heavy impurity 
phase, however, it holds for all n. Therefore the mo- 
ments exhibit multiscaling asymptotic behavior. The de- 
cay coefficients, characterizing the rt-th moment, depend 
on 71 in a nonlinear fashion. This multiscaling behavior 
excludes scaling solutions with sharp tails (stretched ex- 
ponentials decays and faster). While collisions with the 
fluid are sub-dominant, they still lead to corrections to 
the leading asymptotic behavior. 

We turn now to arbitrary spatial dimensions d. The 
impurity velocity changes according to Eq. ( p5| ) with the 
impact direction n chosen randomly. The corresponding 
Lorentz-Boltzmann equation reads 

^Q(v,t)+Q(v,i) = j dnj duQ(u,<) 

x5 V - u + (1 - 9)(u • n) n . (40) 



Moments of the velocity distribution 
M„(t) = / rfvi;"Q(v,i), 



(41) 



can be obtained directly. We focus on the even moments 
of the distribution. Indeed, they satisfy the following 
evolution equation 



dt 



M2n = -{1- {C))M2n 



(42) 



where fi — cos^ 9 — {u ■ n)^, ^ = £_{q, ^) = 1 — (1 — q^)fi, 
and (•) is the shorthand notation for the angular inte- 
gration: (/) = /q I?/i/(/i). Since dn oc sin"^^ 6 d9, the 
(normalized) integration measure is 



B 



1 d- 1 
2' 2 



Vfi = 11 2 (1 — ^) 2' dfi 



(43) 



where B{a, b) is the beta function. For example, (1) = 1, 
and {^) = 1/d. 

From the evolution equations (^), the moments are 
found to decay exponentially with time 



M2n{t) = A'hniO) e 



-{i-{C))t 



(44) 



In particular, 1 — (C) = (1 )/d, and thus, the temper- 
ature decay 8(t) = 8(0) e"^^^'?^-'*/'' is recovered. There 
is a little discrepancy with the exact temperature of 
Eq. ( |l^ ) which depends on the initial fluid temperature, 
Tq. This is a remnant of a transient regime where the 
two velocity scales are comparable. Generally, the time 
dependence is correct. Note also that in one dimension 
i = q^ and hence Eq. (Q) reduces to our earlier result. In 
the infinite dimension limit, /Lt — > 0, and all moments de- 
cay according to e~*. However, in general, the moments 
exhibit multiscaling asymptotic behavior, and knowl- 
edge of the typical velocity is insufficient to fully char- 
acterize the entire velocity distribution. Indeed, writing 
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Man the exponents a„ = (1 - - (0) 

have a nontrivial spectrum. 

The moments directly give a formal exact solution 
of the Fourier transform of the impurity velocity dis- 
tribution (H). We consider isotropic situations where 
G'(k, t) = G(fc^,t) with k = |k|. Expanding the trans- 
form in powers of and substituting the moment result 
(H) yields 



G(k,i) 



n=0 



(2n)! 



M2„(0)e<«">'. (45) 



While this is an explicit solution, it is not too illuminat- 
ing. First, it is in Fourier space, and second, it involves 
the complicated angular averages (^"). 

Nevertheless, it can be shown that the solution re- 
mains a time dependent combination of properly mod- 
ified replicas of the initial distribution. Indeed, ei- 
ther from Eq. (Esl) or directly from the Fourier trans- 
form equation ^^(k, t) + G(k, t) = (G(k^, t)) , the solu- 
tion can be rewritten in the form 



G(fc^^)=e-* V^G„(fc2). 



(46) 



Here, Go(fc^) is the initial Fourier transform, and the 
"building blocks" G„ are obtained from a recursive pro- 
cedure of angular integration G„+i(fc^) = (G„(fc^^)). 



III. THE IMPROVED INELASTIC MAXWELL 
MODEL 



In the IMM, the rates for fluid-fluid and fluid-impurity 
collisions were identical (and therefore, set to unity for 
convenience). For granular fluids, however, the colli- 
sion rate is proportional to the relative velocity [^-^ . 
Therefore, one can improve the Maxwell model by re- 
placing the actual collision rate with an average col- 
lision rate proportional to the average relative veloc- 
ity. Th e simplest ch oice o f the averag e relative veloc- 
ity is V((vi - V2)2) cx V(Ti -|-T2)/2. For fluid-fluid 
and impurity-fluid collisions we thus obtain ^/T and 
\/{T + 0)/2, respectively. Thus, different collision rates 
multiply the collision integrals in the Boltzmann and 
Lorentz-Boltzmann equations (|^)-(^. In the fluid case, 
this overall prefactor merely affects the time dependence 
of the temperature. As will be shown below, in the im- 
purity case, the above phase transition is suppressed in 
the IIMM, although secondary transitions corresponding 
to higher order moments remain. 

Let us again start with the behavior of the tempera- 
ture. The fluid temperature satisfies 



T 



2p{l-p) 



T 



(47) 



Solving this equation, we recover Haff's cooling law 
T{t) = ro[l + t/to]"^, with To the initial temperature 

andio = d/[p(l-p)To'/'] H. 

The corresponding rate equation for the impurity tem- 
perature is 



d T + e 



e 



T 



(48) 



Since we are primarily interested in the temperature ra- 
tio, S = 0/T, we study this quantity directly. It evolves 
according to 



J_ _rf j l + S 

vTdt y 2 



1-9% , 



(49) 



In the inelastic Maxwell model, the gain and the loss 
terms were comparable, both increasing linearly with S. 
Here, in contrast, the loss term, which grows as S^^^, 
eventually overtakes the gain term that grows only lin- 
early with S. Therefore, S —> c where c is the root of the 
cubic equation 
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2p{l-p) 



(50) 



Consequently, there is only one phase, the light impurity 
phase. We note that the ratio c is independent of the spa- 
tial dimension d. Intuitively, since the impurity collision 
rate relatively increases with the impurity temperature, 
the impurity energy dissipation rate increases, thereby 
limiting the (relative) growth of the impurity tempera- 
ture. 




FIG. 3. The equipartition mass rrieq, given by Eq. (^l[). 
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Generally, there is no equipartition of energy except 
for a particular surface in the space {m^rp,rq). Energy 
equipartition occurs when mc = 1. Using the relations 
©"(0)' ^^'^ ^^'^ equipartition mass 




(51) 



Figure 3 plots moq = moq(rp, r,). 

Qualitatively, our findings in the light impurity phase 
of the IMM extend to the IIMM. For example, both ve- 
locity distributions follow scaling forms and the large- 
velocity tails of both distributions are the same. In the 
one dimensional case, explicit expressions for the im- 
purity scaling function are possible, and as the treat- 
ment follows closely that outlined in the light impurity 
phase, we briefly outline the results. In ID, the impu- 
rity velocity distribution approaches a scaling solution 
Q{v,t) T^i/^Q (wTi/^). The corresponding Fourier 
transform reads G{k,T) = g {jk\T^^^y The only differ- 
ence with the above inelastic Maxwell model i s that the 
collision terms are proportional to /3~^ = + c)/2. 
Consequently, Eq. (pOl) generalizes as follows 



(3p{l~ p)zg'{z) + g{z) - (1 + az) e""- g{qz). (52) 



Seeking a series solution of the form (|2l|), leads to the 
following recursion relations for the coefficients 



^ (l-g)g"-i-/3p(l-p) 

1 — — nf3p{l — p) ' 



(53) 



with Aq = Ai = 1. Again, the velocity distribution is a 
combination of powers of Lorentzians as in Eq. (p4|) : 



(54) 



n=2 



where i?„ are linear combinations of the coefRcients An 's 
given by the same expressions as in Sec. II. In particular, 
the large- velocity tail is generic Q{w) ^ w"^. 

Given the algebraic form of the velocity distributions, 
we examine the asymptotic behavior of moments of the 
velocity distribution. Moments of the fluid and the im- 
purity, Ln and M„, respectively, evolve according to a 
straightforward generalization of Eqs.(p7|)-(^, 

—L„+anLn = J2 [■)pH'^-pT'^LjLn-j, (55) 



3=2 
n-2 



— M„ + 6„M„ ■ - qr-'M,Ln-j. (56) 



Here o„(p) and bn{q) are given by Eq. (|2^) and the col- 
lision counters, 

Tp= f dt' Vt, Tg= [ dt' V(r + B)/2, 

Jo Jo 



play the role of time in Eq. ( pSf ) and (^6|), respectively. 
Since the two temperatures are asymptotically propor- 
tional to each other, Tq — > Tp/j3. The fluid moments 

cx f-2'^"/"2 |T7|. In- 



into Eq. IM. 



decay according to L„ (x e" 
serting the asymptotics Ln cx e 
and performing the same analysis as in the IMM we find 
that when Punip) < bn{q), the impurity moments are 
enslaved to the fluid moments, i.e., Af„ oc L„ asymp- 
totically. Otherwise, sufficiently large impurity moments 
behave differently than the fluid moments, viz. L„ oc 
g-hn-r, ^ ^-26„//3o„^ Although the primary transition 
affecting the second moment does not occur, secondary 
transitions affecting larger moments do occur at a series 
of masses, as in Eq. (p2[). In the IIMM nii = di- 
verges, but other masses remain finite. The transition 
masses m„ are found by solving l3a2n{p) = &2n('?) simul- 
taneously with Eq. (50), and then applying Eq. (Q). For 
example, for completely inelastic collisions (r^ = = 0) 
one finds m2 = 1.65. Qualitatively, these transitions im- 
ply that some velocity statistics of the impurity, specifi- 
cally large moments, are no longer governed by the fluid. 



IV. SUMMARY 



We have studied dynamics of impurities in granular flu- 
ids using the inelastic Maxwell model. In general, there is 
breakdown of energy equipartition as two different tem- 
peratures characterize the impurity and the fluid. We 
analyzed two different models. First, we considered iden- 
tical fluid-fluid and fluid-impurity overall collision rates. 
In this case a phase transition, marked by a dissipation 
dependent critical impurity mass, occurs. Breakdown 
of equipartition is moderate in one phase and extreme 
in the other with the asymptotic temperature ratio di- 
verging. Sufficiently light impurities are governed by the 
fluid, and their scaled velocity distribution has similar 
extremal statistics as has the fluid. In one dimension, 
the scaled impurity velocity distribution is given by a se- 
ries containing powers of Lorentzians of all orders (the 
infinite sum truncates for a few special values of the im- 
purity restitution coefficient). Interestingly there is a se- 
ries of increasing transition masses that corresponds to 
divergence of the ratio of velocity moments of decreasing 
order. These masses are always larger than unity and 
the largest such mass corresponds to the temperature. 
When the impurity mass is larger than this critical mass, 
the impurity and the fluid effectively decouple. Next, we 
considered collision rates that account for the tempera- 
ture difference between the fluid and the impurity. In this 
case, the primary transition corresponding to the temper- 
ature is suppressed but further transitions corresponding 
to divergence of higher moments may still occur. 

Many more quantities should be tractable in the frame- 
work of the Maxwell model. For example, one can study 
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velocity correlations as well as the velocity autocorre- 
lation functions. Furthermore, corrections to the lead- 
ing asymptotic behavior, especially in the heavy impu- 
rity phase, can be systematically evaluated. Addition- 
ally, this model can be applied to particles of different 
sizes by modifying the collision cross sections. 

We reiterate that the inelastic Maxwell model exactly 
describes only a mean-field collision process where the 
collision partners arc random and the impact directions 
are chosen according to a uniform distribution. On the 
other hand, it is an uncontrolled approximation when ap- 
plied to real granular fluids. At best, this model may be 
used to approximate sufficiently low velocity moments. 
Still, these results show that the relative collision rates 
play an important role. Comparison with the hard sphere 
Boltzmann equation, Molecular Dynamics simulations, 
and experiment is needed to gauge the utility of either 
of the two approximate models we treated. In practice, 
the impurity problem involves an additional level of dif- 
ficulty. Gathering meaningful statistics requires many 
replica systems with one impurity or alternatively, the 
vanishing volume fraction limit of a mixture. 

This research was supported by DOE (W-7405-ENG-36), 
NSF(DMR9978902), and ARO (DAAD19-99-1-0173). 
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